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A  TECHNIQUE  FOR  SIMULATING  MIXED  NORMAL  AND  UNIFORM  DISTRIBUTIONS 


Because  of  the  suit! variate  nature  of  personnel  problems,  simulation  of 
manpower  systems  requires  techniques  for  generating  sets  of  variables  with 
prespecifiel  Joint  distribution  functions.  For  example,  Sorenson*-,  while 
studying  the  impact  of  manpower  allocation  policies,  has  used  the  Joint  normal 
distribution  to  simulate  the  enlisted  input.  Techniques  used  to  generate 
these  scores  have  been  described  by  Boldt .  Wiskoff,  una  Fitcbr- and  by  Eillier^h 
For  many  problems,  the  normal  approximation  will  be  quite  adequate. 

For  certain  problems,  however,  a  procedure  which  allows  the  generating 
of  both  uniform  mid  normal  variables  is  desirable,  lor  ex;ui>ple,  when  the 
Armed  Forces  Qualification  Test  (AFQ.T)  score  is  included  witn  Array  Classi¬ 
fication  Battery  (AC3)  scores,  the  uniform  distribution  Is  introduced  since 
the  AFQT,  being  a  percentile  measure ,  is  uniformly  distributed.  Conceivably, 
other  percentile  measures  could  be  introduced  aiso.  Iltmce  the  problem  of 
generating  mixed  normal  and  uniform  distributions  can  have  application  in 
personnel  simulations.  Proceiure3  lor  generating  such  distributions  arc 
presented  below  and  their  derivations  outlined. 


CASE  OF  A  SINGLE  UNIFORM  VARIABLE 

Since  procedures  are  available  for  generating  Joint  normal  distribu¬ 
tions  with  specified  means  and  variance -covariance  matrices,  it  will  be 
useful  to  apply  them.  A  brief  review  of  these  procedures  follows.  Means 
will  not  be  cons id-  red  in  this  paper,  since  they  can  be  easily  incorporated 
with  an  additive  adjustment. 

Suppose  it  is  desired  to  generate  a  set  of  variables  Jointly  normally 
distributed  with  variance -covariance  matrix  S.  This  is  accomplished  by 
ge ne rating  a  set  of  variables  which  are  Jointly  normally  distributed  with 
variance -covariance  matrix  equal  to  the  Identity  matrix  to  get  the  matrix  X. 
Then  find  T  such  that  T'T  =  S.  Then  tee  score  matrix  XT  will  contain  vari¬ 
ables  which  are  Jointly  normally'  distributed  with  variance -covariance  mat: ix  S. 
This  procedure  depends  on  the  fact  that  linear  sums  of  jointly  no  really  dis¬ 
tributed  variables  are  themselves-  Jointly  normally’  distributed. 
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The  procedure  suggests  the  possibility  that  to  mix  normals  and  uniforms, 
with  prescribed  variance-covariance  matrices,  the  normals  and  uniforms  could 
be  generated  independently  and  then  transformed  linearly  to  get  the  pre¬ 
scribed  variance -covariance  matrices.  Such  a  procedure  would  indeed  produce 
variables  with  the  proper  variance-covariance  matrices,  but  the  marginal  dis¬ 
tributions  would  be  neither  normal  nor  uniform  since  linear  composites  of 
mixed  distributions  do  not,  in  general,  retain  the  original  distributions. 
Hence,  such  a  procedure  is  not  acceptable. 

The  method  suggested  here  is  to  generate  normal  variables  and  transform 
them  to  a  distribution  whose  variance -covariance  matrix  is  consistent  with 
the  variance-covariance  matrix  desired  when  certain  of  the  variables  are  then 
transformed  bach  into  percentile  measures.  Consider,  for  example,  two 
variables  which  are  Jointly  normally  distributed  with  correlation  p.  Suppose 
one  of  these  variables  is  then  transformed  into  a  uniform  distribution  to  get 
a  variable  which  correlates  r  with  the  other  normally  distributed  variable. 
There  is  clearly  a  functional  relationship  between  p  and  r.  If  the  value  of 
r  were  known  empirically,  then  one  could  find  p,  generate  two  normally  dis¬ 
tributed  variables  with  correlation  p,  transform  one  to  the  uniform  distribu¬ 
tion  to  get  two  variables,  one  normal  and  one  uniform,  with  the  desired 
correlation  r.  The  functional  relation,  whose  derivation  is  outlined  in 
Appendix  A,  has  been  found  and  is 

r  • p  7T  ~ p  *9772*  w 


Hence,  as  one  might  expect,  the  transformation  from  normal  to  uniform  has 
negligible  effect  on  the  correlation.  However,  knowing  the  column  vector 
of  correlations  of  the  normally  distributed  variables  with  the  uniformly 
distributed  variable,  we  compute  (1/.9772)  times  these  correlations  to  get 
the  desired  correlations  for  normally  distributed  variables.  We  then  generate 
all  Jointly  noimally  distributed  variables  und  transform  the  variable  whose 
correlations  have  been  adjusted  into  a  uniform  distribution.  For  example, 
this  procedure  is  suggested  when  an  AFQT  type  variable  is  used  with  ACB  type 
variables . 


CASE  OF  MORE  THAN  ONE  UNIFORM  VARIABIE 

For  many  uniform  variables,  an  analogous  procedure  is  suggested.  In 
this  cose,  though,  one  also  needs  to  find  the  relationship  between  the 
correlation  coefficients  of  two  percentile  variables,  R  say,  which  would  be 
obtained  if  two  normal  variables  with  correlation  R  were  transformed  into 
uniformly  distributed  variables.  This  relationship,  whose  derivation  is 
outlined  in  Appendix  B,  is  given  below. 


Hence,  for  each  value  of  R,  a  value  of  R  con  be  found  which  can  be  used  in 
the  matrix  S. 


DISCUSSION 

Examination  of  equation  (l)  indicates  that  the  correlation  coefficients 
are  not  greatly  affected  by  transformation  from  normal  to  uniform.  Analysis 
of  equation  (2)  indicates  that  the  maximum  difference  between  K  and  K  is 
approx imately  .027.  Hence,  one  can  place  good  reliance  on  studies  in  which 
the  adjustments  were  not  made.  However,  to  add  an  extra  bit  of  representa¬ 
tional  validity  to  future  simulation  studies,  it  is  suggested  that  the 
corrections  be  applied. 
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APPENDIX  W 
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Let  x  and  y  be  Jointly  normally  distributed.  Then 
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